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I. INTRODUCTION 

A Lifshitz pointEHj] is a critical point at which a disordered 
phase, a spatially homogeneous ordered phase, and a spa- 
tially modulated phase meet. In the case of a ci-dimensional 
system with an n-component order parameter, it is called an 
(to, d, rt)-Lifshitz point (or m-axial Lifshitz point) if a wave- 
vector instability occurs in an TO-dimensional subspace. Such 
multi-phase points are known to occur in a wiety of distinct 
physical, systems, including magnetic onestffl, ferroelectric 
crystals^, charge-transfer saltsjffl liquid, cry stalsjl3 systems 
undergoing structural phase transitions!!!! or trying domain- 
wall instabilities]!^ and the ANNNI model.tJ'ta A survey 
covering pthe work related to them till 1992 has been given 
by SelkeJJ which complements and updates an earlier review 
by Hornreich.El Recently there has also been renewed interest 
in the effects of surfaces on the critical behavior at Lifshitz 
pointsMtZl 

From a general vantage point, critical behavior at Lifshitz 
points is an interesting subject in that it presents clear and sim- 
ple examples of anisotropic scale invariance. Epitomized aim 
by dynamic critical phenomena near thermal equilibriumj!a 
and known to occur as well in other static equilibrium systems 
(e.g., uniaxial dipolar ferromagnets), this kind of invariance 
has gained increasing attention in recent years since it was 
found to be realized in many non-equilibrium systems such as 
driven diffusive systems!!!! and in growth processes.! 2 !! 

Systems at Lifshitz points are good candidates-for,studying 
general aspects of anisotropic scale invariance. ElEHl For one 



thing, the continuum theories representing the universality 
classes of systems with short-range interactions at (to, d, n)- 
Lifshitz points are conceptually simple; second, they involve 
the degeneracy m as a parameter, which can easily be varied 
between 1 and d. A thorough understanding of critical behav- 
ior at such Lifshitz points is clearly very desirable. 

The problem has been studied decades ago hY_means of an 
e expansion about the upper critical dimensionErlErEj 

d*(m) = 4+— , m<8. (1.1) 

Other investigations employed the dimensionality expansion 
about the lower critical dimension jLd* (to) = 2 + ^ for 
n > 3, or the 1/n expansion.BEJJnJ Unfortunatel^-the e- 
expansion results to order e 2 one group of authorsBtJtJ ob- 
tained for the correlation exponents rj[ 2 and 77/4 and the wave- 
vector-exponent j3 q are in conflict with those of Sak and 
GrestrJ for the cases to = 2 and to = 6. 

In order to resolve_this long-standing controversy, Mer- 
gulhao and CarneircoEj recently presented a reanalysis of 
the problem based on renormalized field theory and dimen- 
sional regularization. Exploiting the form of the resulting 
renormalization-group equations, they were able to derive var- 
ious (previously given) general scaling laws one expects to 
hold according to the phenomenological theory of scaling. 
However, their calculation of critical exponents was limited 
in a twofold fashion: They treated merely the special cases 
to = 2 and to = 6, in which considerable simplifications oc- 
cur. Their results for 77^ and 77/4 to order s2_ agree with Sak 
and Grest'sElbut disagree with Mukamel'sS Second, the ex- 
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ponent f3 q (an independent exponent that does not follow from 
these correlation exponents via a scaling law) was not consid- 
ered at all by them. Thus it is an open question whether Sak 
and Grest's or Mukamel's 0(e 2 ) results for j3 q with m = 2 
and to == 6 are correct. Furthermore, for other values of m, 
the published 0{( 2 ) resultscjo for the exponents 7^2, rju, and 
(3 q remain unchecked. 

It is the purpose of this work to fill these gaps and to deter- 
mine the e expansion of the critical exponents 77/2, rju, and (3 q 
for general values of to to order e 2 . 

Technically, we shall employ dimensional regularization in 
conjunction with minimal subtraction of poles in e. This way 
of fixing the counterterms appears to us somewhat more con- 
venient than the use of normalization conditions (as was done 
in Refs. M and 30). In order to overcome the rather demand- 



ing technical challenges, we have found it useful to work di- 
rectly in position space. Thus the Laurent expansion of the 
distributions to which the Feynman graphs of the primitively 
divergent vertex functions correspond in position space must 
be determined to the required order in e. 

The source of the technical difficulties is that these Feyn- 
man graphs, at criticality, involve a free propagator G(x) 
which is a generalized homogeneous rather than a homoge- 
neous function, because of the anisotropic scale invariance of 
the free theory. While such a situation is encountered also in 
other cases of anisotropic scale invariance, the scaling func- 
tion associated with G(x) turns out to be a particularly com- 
plicated function in the present case of a general (to, d, n)- 
Lifshitz point. (For general values of to, it is a sum of two 
generalized hypergeometric functions.) 

In the next section, we recall the familiar continuum model 
describing the critical behavior at a (to, d, n)-Lifshitz point 
and discuss its renormalization. In Sec. |l| details of our cal- 
culation are presented, and our results for the renormalization 
factors are derived. Then renormalization-group equations are 
given in Sec. £v|, which are utilized to deduce the general scal- 
ing form of the correlation functions, to identify the critical 
exponents, and to derive their scaling laws as well as the an- 
ticipated multi-scale-factor universality. This is followed by 
a presentation of our e-expansion results for the critical ex- 
ponents rii2, 7]4i, and (3 q . Sec. contains a brief summary 
and concluding remarks. Finally, there are two appendices to 
which some computational details have been relegated. 



II. THE MODEL AND ITS RENORMALIZATION 

We consider the standard continuum model representing the 
universality class of a (m, d, n)-Lifshitz point with the Hamil- 
tonian 

H{4>] = \ J d d x[ Po (V|| <j>) 2 + a (A||0) 2 

+ (V ± ^) 2 +r o 2 + ^|0| 4 } . (2.1) 

Here 4>{x) = (<pi, . . . ,<f> n ) is an n-component order- 
parameter field. The coordinate x 6 M. d has an Tri- 



dimensional parallel component, xn, and a (d — Tri- 
dimensional perpendicular one, x±. Likewise, V|| and Vi de- 
note the associated parallel and perpendicular components of 
the gradient operator V, while A|| means the Laplacian Vy 2 . 
At the level of Landau theory, the Lifshitz point is located at 

p = T = 0. 

The Hamiltonian is invariant under the transformation 



00 



a 4 (jQ , 
a m u Q 



Po 



■ x± , <j> -> 
a 2 Po , t 



-m/2 



TO 







(2.2) 



Thus, appropriate invariant interaction constants are 
wooo - ™/ 4 , poco -1 / 2 , and to, and the dependence on the 
parallel coordinates is through the invariant combination 

-1/4 
(T ' X\\- 

Dimensional analysis yields the dimensions [.]: 



M=M 2 , bo] = [ct ] 1/2 /J > 
N = h)] W V with e = d*(m) - d 
i(x)j = 



where p is an arbitrary momentum scale. 
Let 



Lr ii,...,i N \. x li 



,x N 



h (x 1 )...(j) lN {x N )} c 



(2.3) 



(2.4) 



denote the connected iV-point correlation functions (cumu- 

lants) and i (xi,...,xn) the corresponding vertex 

functions. Using power counting one concludes that the ul- 
traviolet (uv) singularities of these functions can be absorbed 
through the reparameterizations 
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where 



A 



d,m 



$d—m $m — 



4 7^/2 



r(±^a)r(m/2) 



(2.6) 



is a convenient normalization factor we absorb in the renor- 
malized coupling constant. Here 



S D = 



2tt j / 2 
T(D/2) 



(2.7) 



is the surface area of a D-dimensional unit sphere. 

The quantities to c and poc correspond to shifts of the Lif- 
shitz point. In our perturbative approach based on dimen- 
sional regularization they vanish. If we wanted to regularize 
the uv singularities via a cutoff A (restricting the integrations 
over parallel and perpendicular momenta by |gn| < A and 
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\q±\ < A), they would be needed to absorb uv singularities 
quadratic and linear in A, respectively. 

In the renormalization scheme we use, the renormalization 
factors Ztj,, Z a , Z T , Z p , and Z u , for given values of the pa- 
rameters e, n, and m, depend just on the dimensionless renor- 
malized coupling constant u; that is, they are independent of 
a, t, and p. This follows from the fact that the primitive di- 
vergences of the momentum-space vertex functions r( 2 )(q) 
and f^ 4 ) (q x , . . . , q x ), at any order of uq ero~ m / 4 , are poles in 
e whose residua depend linearly on q± 2 , po q\\ 2 , oo q\\ A , and 
t in the case of the former and are independent of these mo- 
menta and mass parameters in the case of the latter. Subtract- 
ing these poles minimally as usual implies that these renor- 
malization factors differ from 1 through Laurent series in e: 



z u = i + ^2 a S ( u ; m > n ) e ~ 



P =i 

oo r 



r—1 p—1 



u, t, a, p. 



(2.8) 



III. OUTLINE OF COMPUTATION AND PERTURBATIVE 
RESULTS 

We shall compute the leading nontrivial contributions to 
these renormalization factors. In the cases of Z$, Z a , and 
Z p , whose 0(u) contributions vanish, these are of order u 2 ; 
for Z u and Z T they are of first order in u. 

To this end we expand about the Lifshitz point, using the 
free propagator 



G(x) 



J-{*iw x \\ +q±-x±) 

cto q\ + q± 



(3.1) 



Here the (dimensionally regularized) momentum-space inte- 
gral is defined through 



d d ~ r 



q 



d m q 

(2ir) m J (2ix) d - m 



(3.2) 



q« q± 



Let ?'|| = |a;||| and r± = \x±\. Then the free propagator can 
be written in the scaling form 



G( ; r)=r ± - 2+£ <i>( ( 7o- 1/4 r|| r^ 1 ' 2 ) 



with 



$(u) = to, d) 



qt + qi 



(3.3) 



(3.4) 



where v 6 R m is a vector of length v and arbitrary orienta- 
tion, while e± means the unit vector x±/r±. Note that the 
scaling function $ depends parametrically on to and d. For 



the sake of brevity, we will usually suppress these variables, 
writing Q(v; to, d) only when special values of to and d are 
chosen or when we wish to stress the dependence on these 
parameters. 



The integration over q ± in (3.4) yields 



= (2tt)- 



q^-^^-MO 



(3.5) 



Upon introducing spherical coordinates q = \q«\ and 
ft( m ) = (0i,...,0 m _i) for qp with dn^ = 
sin" 1-2 9 m -i d6 m - \ df^™ 1 " 1 ), one can perform the angular 
integrations. This gives 



$(i>) 



v 2 
(2 7r) d / 2 



dqq 



2-e 



:MKtE_!(9 2 ) • (3.6) 



The integral remaining in (3.6) can be expressed as a combi- 
nation of generalized hypergeometric functions 1F2 (see Ap- 
pendix |a|). For special values of to and d, the result reduces 
to simple expressions, which we have gathered in Appendix 

The leading loop corrections to the vertex functions T'- ' 
and r^ 2 ) at the Lifshitz point are given in position space by 
the graphs x <(^^P y and x <^^> y , which are propor- 
tional to G 2 (x — y) and G 3 (x — y), respectively. Hence we 
must determine the Laurent expansion of these distributions. 
To this end we set = 1 and consider the action (G s , tp) 
of G s (x) for s = 2,3 on a test function f(x). We substi- 
tute ( |3.3| ) for G and use spherical coordinates (riufiii) and 
(rj_ , fijj for the parallel and perpendicular components of x, 
writing ip(x) = <p(rt\ , fin ; r±, fl±). We thus obtain 

(G S ,V) =(r ± - s ^&(i 1 r ± - 1 / 2 ),^ 



J d d xr ± - s ^ 2 -^ 



-1/2) 



tp(x) 



^ i r i - > N+!^( Ii ) ] ( 3.7) 



where the functions ip s (x±) = ^p s (r±,il±) are defined 
through 

i/>.{x ± ) = /d m X||^(r||)^(r|| ^I,n f ,r ± ,n ± ) . (3.8) 



The final result in (3.7) is the linear functional 

(Y^ - ^ 2 " 6 )" 1 ^ , t/> s ). Generalized functions such as rj_^"'' 

and their Laurent expansions are discussed in Ref. |T} Let 

ip(x±) = be a smooth (C°°) test function on 

R d- m and 

^(rijE-^ frfn^fri^i) (3.9) 
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its spherical average. Then we have 



-s(2-e)-M 



S d - m /drr 3 - 3 '^- 1 )^^) 
o 



3-2 S +e( S -l) T^J 



(3.10) 



Here r+ A is a standard generalized function in the notation of 
Ref. glj. Its Laurent expansion about the pole at A = — p = 
— 1, —2, . . . reads 



(_l)p-i 5(p-i)( r ) 



0-1)! A + p 
where the generalized function r+ ~ p is defined by 

p-2 



0(A+p), (3.11) 



(r + -",v»(r)) - ydrr + -f 



(P-1) 



^w-E-^ (J) (o) 

3=0 3 ' 



(3.12) 



Using these results, the leading terms of the Laurent expan- 
sions of (G s , p) can be determined in a straightforward man- 
ner. However , it s hould be noted that the functions ips(x±) 
introduced in (3.8 ) are not a priori guaranteed to have the usu- 



ally required strong properties of test functions (continuous 
partial derivatives of all orders and sufficiently fast decay as 
— * oo). In particular, one may wonder whether the de- 
pendence on the variable rn \frZ of y in (3.8) does not imply 
that derivatives such as W±tp s become singular at the origin. 
Closer inspection reveals that this is not the case since the 
problematic term ~ r^ 1 involves the vanishing angular in- 
tegral Jdf2|| xii (p(. . .). 
One obtains 



(3.13) 



and 



Sd—m 



V>3 n ""(0) 



4e 



(3.14) 



From its definition in (3.8) we see that the residuum ^2(0) on 
the right-hand side of ( 3.13 ) reduces to a simple expression 
oc <ys(0). We thus arrive at the expansion 



9p*± = JoArn,d*) §{x)+0{e0) 
where J 0:2 is a particular one of the integrals 



(3.15) 



J p . s {m,d)= Jv m - 1+p $ s {v-m,d)dv 





(3.16) 



In order to convert the Laurent expansion ( 3.14 ) into one 

for G 3 (x), we must compute xps (0). This in turn requires 
the calculation of the following angular average: 



o^n r \\v r > n \\> r > n ±) 



r=0 



-Hy ,Q 



-(^n-Vii) p(0) + (e x • V ± r V(0) 



4! 



r|| 4 (A||V)(0) | (A ± p)(0) 
4m(m + 2) d — m 



■O11 



(3.17) 



Using this in conjunction with ( 3.14 ) gives 

G 3 {x) _ J 4 ,3Kd*)A|| 2 (5(a;) 



16 m(m + 2) e 
J , 3 (m,d*)A_ L 6(x) 
4 (d* —m)e 



O(e ) 



(3.18) 



In order to compute the 0(u 2 ) term of Z p , we consider 
the two-point vertex function with an insertion of the operator 
i J d d x(S7 \\4>) 2 (to which po couples). We represent such an 
insertion by the vertex -h- . At the Lifshitz point r = p = 0, 
the leading nontrivial contribution to this vertex function is 
given by the two-loop graph x . The upper line in- 

volves the convolution 



- (V||G* VyG)(x) = 
where 

= S(f ; m, d) = 



g..2 e i(g||-u+9±-ej.) 



('HI 



(3.19) 



(3.20) 



is the analog of the scaling function (cf. (3.3)). Proceed- 
ing as in the case of the latter, one obtains 



E(v) 



1 



2(271-)* 



2(2 7r) d / 2 



„ d—m—2 \r 



: (<7||V"' 911 



dq q 2 - £ 3^2(qv) K d- m -4 (q z ) 



(3.21) 



The remaining single integral can again be expressed in terms 
of generalized hypergeometric functions. The corresponding 
general expression, as well as the simpler ones to which this 
reduces for special values of m and d, may be found in Ap- 
pendix ^). 
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The required graph x is proportional to the dis- 

tribution 



D(x) = -G 2 {x) (ViiG* Vi|G)(aj) 



(3.22) 



whose pole term can be worked out in a straightforward fash- 
ion by the techniques employed above. One finds 



G 2 {x) (ViiG* V|iG)(x) Ji(m, d*) Am*(x) 



A 



4m 6 



with 



•0(e°) 
(3.23) 

(3.24) 



A convenient way of computing the renormalization factor 
Z T is to consider the vertex function I^ 2 ' 1 ) with a single inser- 

V 

tion of the operator | <fi(y) 2 , which we depict as -i- . Its one- 
y, 

loop contribution C^) is proportional to G 2 (x — y). Hence 
x 

its Laurent expansion follows from that of the latter quantity. 

Let us introduce coefficients b L (m) for the leading non- 
trivial contributions to the renormalization factors Z L , writing 
these in the form 



Z u = l + b u (m)^-+0(u 2 ) 



, . . n+2 u 9 
1 + b T m — - - + (u 2 ) , 
3 e 



(3.25) 



(3.26) 



and 



„ , , . n+2 u 2 

Z q = 1 + b, m — - + O (u 3 ) 
3 e 



<; = <f>,a,p. (3.27) 



From the pole terms of G 2 (x — y) given in (3.15) one easily 
deduces that 



b u (m) = 3 b T (m) = - J 0i2 (to, <f ) 



(3.28) 



The pole terms proportional to A±5(x), A\ 2 5(x), and 
A || S(x) of the two-loop graphs considered above are ab- 
sorbed by counterterms involving the renormalization factors 
Zf, Z a = Z fj Z$, and Z p = Z Z& Z a 1/ 2 res pectively. Uti- 
lizing the Laurent expansions (3.18 ) and (3.22), one finds that 
their coefficients are given by 



b${m) 



1 



1 J ^(m,d*) 



24 d* — 771 A n 



b a {m) = — 



J4,3(m, d*) 



96 m(m+2) A d% 



(3.29) 



(3.30) 



and 



1 h(m,d*) 



8m A d . 



The coefficients b a and b p are related to these via 



and 



b a (m) = 6 CT (m) - b^(m) 



b p (m) = b p (m) - - b<f,(m) - ^ b a (m) 



(3.31) 



(3.32) 



(3.33) 



IV. RENORMALIZATION GROUP EQUATIONS AND 
e-EXPANSION RESULTS 



The reparameterizations (2.5) yield the following relations 
between bare and renormalized correlation and vertex func- 
tions 



G^(x h x ± ) = Z^ 2 G^(x h x ± ), 

rW(x h x ± ) = z 4> - N / 2 rW(x h x 1 _) 



(4.1a) 



(4.1b) 



where xn and x± stand for the set of all parallel and perpen- 
dicular coordinates on which G^- 1 andr( w ^ depend. For con- 
ciseness, we have suppressed the tensorial indices 
of these functions and will generally do so below. 

Upon exploiting the invariance of the bare functions under 
a change u — ► p(£) — p I of the momentum scale in the usual 
fashion, one arrives at the renormalization-group equations 



' N 



_ N 



G w = 

w ren J 

r W = 



(4.2) 



(4.3) 



with 



P d p + U d u - 7]o- a d a 

-{2 + 7] T )Td T - (1 +7] p )pd p , 



where the beta and eta functions are defined by 

(i u = jj,dp\ u= -u [e + ri u (u)} 

and 

r) t = p dp\ Q In Z, , l = <p, a, p, r, u , 



(4.4) 



(4.5) 



(4.6) 



respectively. Here <9 M | Q means a derivative at fixed bare vari- 
ables po, po, (To, and to. Owing to our use of the minimal 
subtraction procedure, the functions r) L can be expressed in 
terms of the residua a t) i [u\ m, n) as 



da^i 
du ' 



a, p, t, u. 



(4.7) 
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To solve the renormalization-group equations (4.1) via is positive. 



characteristics, we introduce flowing variables through 

1^,U{1)=P U [U(1)\, U(\) = U, 



The solutions to the other flow equations, (4.9)-( 4.1 1 ), can 
be conveniently written in terms of the anomalous dimensions 
(4.8) ?/* = t] l (u*) and the renormalization-group-trajectory inte- 
grals 



and 



=-[l + VpW\P, P(V=P, 



e-nT(Z) = -[ 2 + vAu)]T, f(l) = r 



(4.9) 



(4.10) 



(4.11) 



The flow equation (4.8) for the running coupling constant 
u(£) can be solved for I to obtain 



E L [u, u] = exp < 



u(l) 



d,r 



which approach nonuniversal constants 

E *( u ) = E L {u*,u) , L=<j),a,p,T. 

in the infrared limit I — ► 0. 
One finds 



a, p,r. 



(4.15) 



(4.16) 



ln£ = 



dx 

0u{x) 



(4.12) 



For e > 0, the beta function /3 u (u) is known to have a non- 
trivial zero u*, corresponding to an infrared-stable fixed point. 
Expanding about this fixed point gives the familiar asymptotic 
form 



u(£) = u* + (u-u*)F* + 0(£ 2u ") 
in the infrared limit £ — ► 0, where 

W« = -3— (tt ) 
du 



(4.13) 



(4.14) 



cf(^) = r< E a [u(£),u]a 
£'< E;(u) a, 



and 



p{£) = 



--(£) = £-^ +r >^ E T [u{l), 



Solving the RG equation ( 4.1a ) in terms of characteristics yields 

-1 N/2 



G^(xpx±;p,T,u,a,p) = 

N/2 

CT™/ 4 E,f,{u,u) 
To obtain the second equality, we have used the relation 



Ft 



E<f,(u,u) 



G^(x\\,xx;p,T,u,a,p,£) 



/i 2 ^ 2 /cr) 1 ^ 4 xu, p,£x±; p, T,U, 1, 1 



GW(x||,»_l;p,t,'u,ct,/x) = 



AT/2 



Gffl (<t" 1/4 /2 1/2 x|| , /2 xx; p, f, 72, 1, 1) 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



(4.21) 



implied by our dimensional considerations (2.3). 

Let us assume that the function gJw on the right-hand side of ( 4.2C ) has a nonvanishing limit u — > u* for e > 0. This 
assumption is in conformity with, and can be checked by, RG-improved perturbation theory. We choose £ = £ T such that 



f(£ T ) = ±1 for ± t > 



(4.22) 



and consider the limit r — * 0±. To write the resulting asymptotic form of gJw in a compact fashion, we introduce the 
correlation-length exponents 



Vl2 = 



1 

2 + »£ 



(4.23) 
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and 



2 + 



the crossover exponent 



as well as the correlation lengths 



and 



4(2 + V * T ) ' 
ip = v n (I+77*) , 

i^ = vr x i T ^vr x [e;{u)\t\]-^ 



jj? i 2 



1/4 



p- 1 ' 2 [E;(u)a} 1/4 [E;(u)\r\Y 



In terms of these quantities the asymptotic critical behavior of Grer? becomes 



E* 



N/2 



Gi N) 



X|| x 



with 



(4.24) 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



The result is the scaling form expected according to the 
phenomenological theory of scaling. As it shows, the scal- 
ing function Q±' is universal, up to a redefinition of the 
nonuniversal metric factors associated with the relevant scal- 
ing fields, i.e., E*, E*, E*, and E^. (Note that E%, whose 

change would affect the overall amplitude of G± , as usual 
corresponds to a metric factor associated with the magnetic 
scaling field; see, e.g., Ref. p2| .) 

The correlation exponents rji2 and 7^4 are given by 



and 



V12 = % 



2 + r?* 



(4.30) 



(4.31) 



This can be s een e ither by taking the Fourier transform of the 
above result (4.28 ) with N = 2 or else by solving directly the 
renormalization-group equation of Trcn^n ,Q±)- In order to 
identify the wave-vector exponent f3 q , we utilize the scaling 
form 



mi(q\\,q±; T ,p, u ) ~ M 7 ^±(q\\^\\,q±^±;p\T 



(4.32) 



of the inverse susceptibility f ^ and argue as in Ref. |23|: On 
the helical branch T^\(p) of the critical line, the inverse sus- 
ceptibility vanishes at q c = (qjj, 0) 7^ 0. Hence in the scal- 
ing regime, the line Th e i(p) is determined by the zeroes of the 



scaling function T(p, 0, g). Denoting these as p c and g c , we 
obtain the relations 



q\\=Pc£\\ 1 ~Pc\t\ 



and 



which yield 



with 



P = 0c\T\ 



111 



P q = — = 



2 + 



ip 4(1 + 77;) 



(4.33) 



(4.34) 



(4.35) 



(4.36) 



where the last equality follows upon substitution of (4.25) and 
(4.24) for <p and ^4, respectively. 



To com pute the exponent functions (4.6) and the beta func 



tion (4.5b, we insert the residua of the renormalization factors 



(3.25)-(3.27) into (4.7) and express b T in terms of b u using 



(3.28). We thus obtain 



77+2 

% ( u ) = -2 — b,(m)u 2 + 0(u 3 ) 7 <; = </>,*, p, (4.37) 



rj T {u) = - i b u (m) u + 0(u 2 ) 



(4.38) 



and 



7 



n+8 
~9~ 



b u (m)u + 0(u 2 ) 



(4.39) 



From the last equation we can read off the e expansion of 
u*, the nontrivial zero of (3 U : 

9 



n+8 b u (m) 



+ 0(e 2 ) 



(4.40) 



Evaluation of the above exponent functions at this fixed- 
point value gives us the e expansions of the anomal ous d i- 
mension s r\* . Substituting these into the expressions (4.23)- 
(|4.25|), (OS), (Oil), and (OJ) for the critical exponents 



yields 



1 



n+2 



4(n+8) 
27(n+2) & (to)-6 ct (to 



0(e 2 ) 



(4.41) 



(n+8) 2 2b u (mY 



0(e 3 ) 



0(e 3 



27(n+2) e 2 
(n+8) 2 T 



r 0.02152 for m = 1, 

0.02195 for m = 2, 

0.02231 for m = 3, 

0.02263 for m = 4, 

0.02290 for m = 5, 

I 0.02313 for to = 6, 
(4.42) 



«; 2 = -2 



27(n+2) 6 (m) 



(n+8) 2 6„( m ) 



0(e 3 ) 



27(n+2) 



n+£ 



e 2 < 



«;4 = -4 



27(n+2) b a (m) 



(n+8) 2 6 u (m) 2 



0(e 3 ) - 



27(n+2) 
(n+8) 2 



e 2 < 



e 2 +G(e 3 ) 






' 0.01739 


for to = 


1. 


0.01646 


for to = 


2, 


0.01564 


for m 


3, 


' 0.01488 


for m 


4, 


0.01418 


for m 


5, 


. 0.01353 


for to = 


6, 






(4.43) 


e 2 +0(e 3 ) 






' 0.00827 


for to 


1, 


0.01097 


for m 


2, 


0.01334 


for to = 


3, 


' 0.01548 


for to = 


4, 


0.01743 


for to = 


5, 


, 0.01920 


for to 


6, 



VII 



27(n+2) b$(m)-2b p (m)+b a (m) ^ 2 



(n+8) 2 



(4.44) 



0(e 3 ) 



1 + C(e 3 



27(n+2) 
(n+8) 2 



6 u (m) 

f 0.02781 for to =1, 

0.05487 for to = 2, 

0.07856 for m = 3, 

0.09980 for to = 4, 

0.11904 for to = 5, 

1 0.13658 for to = 6, 

(4.45) 



e 2 < 



and 



ft 



27(n+2) 6 p (to) 

^2 



-b a (m) 



(n+8) 2 



G(e 3 



b u (m) 



0(e 3 ) 



27(n+2) 
(n+8) 2 



r 0.00852 


for 


TO 


= 1. 


0.02195 


for 


TO 


= 2. 


0.03370 


for 


m 


= 3. 


0.04424 


for 


m 


= 4. 


0.05379 


for 


m 


= 5. 


^ 0.06251 


for 


m 


= 6. 








(4.46) 



We have expressed the 



cie nts b u ( to), b , 



b a (m 



results 
and 



in terms of the coeffi- 
b p (m), which according 



to ( 3.28| )-(3.31 ) are proportional to the integrals Jo,2(ni, d* 
Jo ; s(m,d*), J4,3(m, d*), and h (m, d *), re spectively. These 
integrals are defined by ( 3.16 ) and ( 3.24 ). The first one 
of them — the one-loop integral Jo,2(m, d) — is analytically 
computableLJ for general values of d and to. The result is 

O-2-e 



Jo,2(w, d) 



(2k) 



4 



giving 



bu(m) = - 



r(2-f- £ ) 

T(2-e) 



r(f) 



(27T) 



4+^ 



(4.47) 



(4.48) 



The fixed-point v alue that results when this value of b u (m) 
is inserted into ( 4.40| ) is consistent with the one found in 
calculations based on Wilson's momentum-shell integration 
mefhod.Ej 

The integrals Jo,s(m, d*), J^^m, d*), and Ii(m, d*), and 
hence the coefficients 6^(to), b a (m), and b p (m), can be cal- 
culated numerically for any desired value of m, using the ex- 
plicit expressions for the scaling functions $(u;m, d*) and 
H(v; to, gT) given in (A4) and (A5) of Appendix 0. (As dis- 
cussed there, the numerical evaluation of these integrals for 
general values of to requires some care because Q(v; m, d*) 
is a difference of two terms, each of which grows exponen- 
tially as v — * oo.) In this manner one arrives at the values of 
the e 2 terms given in the second lines of ( 4.42 )-( 4.46). 



In Fig. |T] the coefficients of the e 2 terms of some of these 
exponents are depicted for the scalar case, n = 1. As one 
sees, they have a smooth and relatively weak rn-dependence, 
especially for r}i% and r/u. 



0.06 



.04 



. 02 



-0 . 02 
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FIG. 1. Coefficients of e 2 terms of the exponents r\i2 (triangles), 
■qu (stars), and j3 q (squares) for n = 1. 

In the special cases m = 2 and m = 6, the functions 
$(t> ; m, d* ) and S(t>; m, c?*) become sufficiently simple [see 
(A6)-( A8|], so that the required integrations can be done an- 
alytically. This leads to 



M2) 



M2) 



1 1 

54 (4tt) 



1 1 

162 (4tt) 8 



Sp(2) = 



1 1 

18 (4tt) 8 



M6) 



6.(6) 



16 l-31n| 
~9 (4tt) 12 



14 1 

81 (47T) 1 



and 



M6) = 7T 



1 



61n| 



9 (4tt) 15 



(4.49a) 



(4.49b) 



(4.49c) 



(4.50a) 



(4.50b) 



(4.50c) 



If these analytical ex press io ns for the co effici ents are inserted 
into the expansions ( [4.43 ), (4.44), and (4.46) °ffH2, T]i4,, and 
q with m = 2 and to = 6, then Sak and Grest'sHa results for 
those two values of to are recovered (which in turn agree with 
Mergulhao and Carneiro'sEj findings for r\i2 and 77/4). 

As was mentioned already in the Introduction, these results 
for m — 2 and to — 6 disa gree wi th M ukameFsEj. More 
generally, our £>(e 2 ) results ( 4-.42|)— ( p~46b , for all values of 
m = 1, . . . , 6, turn out to be at variance with the latter au- 
thor's,— The case to = 1 was also studied by Hornreich and 
BruceE.3, who calculated 77(4(771=1) and f3 q (m=\) to order e 2 . 
Their results agree with Mukamel's and hence diagree with 
our's. 



Upon extrapolating the series expansions ( 4.42 )-( 4.46 ) one 
can obtain exponents estimates for three-dimensional systems. 
Unfortunately, there exist in the literature only very few pre- 
dictions of exponent values produced by other means with 
which we can compare our's.cl Utilizing high-temperature 
series techniques, Redner and Stanleys found the estimate 
(3 q = 0.5 ±0.15 for the case of a uniaxial (to, d, n) = (1, 3, 1) 
Lifshitz point. This is in conformity with the value j3 q ~ 
0.519 on e gets by setting e = 1.5 in the corresponding m=l 
result of ( 4.46 ). A mare recent high-temperature series analy- 
sis by Mo and FererEa yielded 2 (3 q ~ 1. For the susceptibility 
exponent 



the correlation exponent 144, and the specific-heat exponent 



at 



2 — m via — (d — tri) vn 



(4.52) 



of the (m, d, n)=(l, 3, 1) Lifshitz point these authors ob- 
tained the results 7i = 1.62 ± 0.12, 4i/ M = 1.63 ± 0.10, and 
ui = 0.20 ± 0.15. Utilizing these numbers to compute 77/4 
via the scaling law implied by ( 4.5 1 ), ^4 = 4 — 72 /V;4, yields 
r]i4 ~ 0.02 ± 0.5. This may be be compared with the value 
~ —0.019 one finds from (4.44) upon setting e = 1.5. r— ■ 
As a further quantity for which Mo and Ferer's resultst3 
yield an estimate that can be compared with our 0(e 2 ) results 
we consider the ratio Pi/ji. Substituting their exponent val- 
ues into Pi = (2 - ai - 7;)/2 yields # = 0.09 ± 0.135 
and Pi/ji — 0.055^q;q94. From the asymptotic form (3.3) of 

Crcn _ ^ ^ one reads off the scaling law 



Pi = — {d - m - 2 + r]i2j + — m 



(4.53) 



which may be combined with relation (4.51) for 7; to conclude 
that 



7/ 

We now set m 
lives vu/ v n 



2(2 - m2 ) 



(4.54) 



n = 1 and e 
0.488 and r/ ;2 



1.5 in 04.420 and ( |4.43| ). This 
0.039. 



Then we insert these 
numbers into ( J4.54[ ) with d = 3, obtaining /3;/7/ ~ 0.134. 

There also exist Monte Carlo estimates of-exaonents for the 
case of a (to, d, n)=(l, 3, 1) Lifshitz point£jE3 The more re- 
cent ones_& = 0.19±0.02and7 ; = 1.4± 0.06, due to Kaski 
and SelkeEZI, give /3 ; / 7/ = 0.136 ± 0.02. In view of the fact 
that the importance of anisotropic scaling and its implications 
for finite-size-effects in systems exhibiting anisotropic scale 



invariancel£3i±J has been realized only more recently, it is not 
clear to us how reliable these Monte Carlo estimates may be 
expected to be. Note, on the other h and, that the coefficients 
of the e 2 terms of the series ( 4.42 )— ( J4.46 ) are all truly small. 
Thus it is not unlikely that the values one gets for d = 3 by 
naive evaluation of these truncated series are fairly precise, at 
least for to = 1. (The e 2 -corrections of these exponents grow 
with to because of the factor (d* - 3) 2 = (1 + y) 2 -) 



li = v n (2 - rji 2 ) = 144(4 - 77; 4 ) 



(4.51) 



V. CONCLUDING REMARKS 

We have studied the critical behavior of d-dimensional sys- 
tems at m-axial Lifshitz points by means of an e expansion 
about the upper critical dimension d* = 4 + y . Using mod- 
ern field-theory techniques, we have been able to compute the 
correlation exponents r)i2 and 77/4, the wave-vector exponent 
(3 q , and exponents related to these via scali ng law s to ord er e 2 . 
The resulting series expansion^- given in (4.42 )-( 4.46| ), ct 
rect earlier results by MukamelEj and Hornreich and Bruce 
for the special values to = 2 and m = 6, we recovered Sak 
and Grest'sc3 findings. 

To clarify this long-standing controversy, it proved use- 
ful to work directly in position space and to compute the 
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Laurent expansion of the dimensionally regularized distribu- 
tions associated with the Feynman diagrams. There are two 
other classes of difficult problems where this technique has 
demonstrated its potential: field theories of polymerized (teth- 
ered) membxapesc3c3 and critical behavior in systems with 
boundaries. l3cj In the present study an additional complica- 
tion had to be mastered: The free propagator at the Lifshitz 
point, which because of anisotropic scale invariance is a gen- 
eralized homogeneous function rather than a simple power of 
the distance \x — x'\, involves a complicated scaling func- 
tion. For powers and products of simple homogeneous func- 
tions, a lot of-mathematical knowledge on Laurent expansions 
is available. EJ Unfortunately, the amount of explicit mathe- 
matical results on Laurent expansions of powers and products 
of generalized homogeneous functions appears to be rather 
scarce. Since we had no such general mathematical results at 
our disposal, we had to work out the Laurent expansions of 
the required distributions by our own tools. 

Difficulties of the kind we were faced with in the present 
work may be encountered also in studies of other types of sys- 
tems with anisotropic scale invariance. Hence the techniques 
utilized above should be equally useful for field-theory analy- 
ses of such problems. 
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APPENDIX A: THE SCALING FUNCTIONS AND H(w) 



The scaling func tions and 5( V) i ntroduced respec- 

tively throug h d3.3| )-( [3~4|) an d (3. 19 )— ( 3.2C| ) are given by sin- 
gle integrals (|3.6|) and (]3.21|) of the form 



i(v) = v 



dq q 2 



e J M (gu)IW) 



(Al) 



This is a standard integral]?! which for arbitrary values of 
its parameters /i and v, can be expressed in terms of gener- 
alized hypergeometric functions 2^3- For the special values 
ji=f-landi/ = l- f- 
is needed, it simplifies, giving 



I or v = — ?f — I for which it 



to, d) = — e+7^r, 

2 2+m 7T 4 



and 



S(u; to, d) 



23+m n 4 



r(i-f) 
r(^) 

£(¥) 
r (f) 



2' 2' 4 '64/ 4 r(l+f ) V 2 2 '2' 4 ' 64 



1F2 



r(|-|) F (3 e.3 my 
1^2 77 — 77, r> J-+-T) 777 



1-e 1 m_ u 4 \ v 2 r(l-f) 

~ ' 2' T' 64 J ~ Tr72±™y 



e 3 2+m v 4 

' il 2 : 2' i - : (iT 



(A2) 



(A3) 



At the upper critical dimension, i.e., for e = 0, this becomes 

1 



$(u;m,4+^) = 



2 5+m 7T" 



1 2+771 V \ 3rn _ 9 m (V \ 



(A4) 



and 



Hfu; to, ' 



. 6+TT 

2 6+m 7T~ 



2— v^Fu 2 "- — ,^ r iF 2 1; -, 



3 2+to f 



4 ' 64 



(A5) 



respectively, where the I„(.) are modified Bessel functions of the first kind. 



In the special cases m = 2 and m = 6, these expressions 
reduce to simple elementary functions: One has 



S(u; 2, 5) = -$(«; 2, 5), 



(A7) 



*(t>;2,5) = 



(4tt) 2 



(A6) 



3(u;6,7) = 



l-(l + x)e" 



(2tt) 3 tj 4 



(A8) 



10 



and 



H(u;6,7) 



(An) 



1 1 / 

7r) 3 V 2 V 



(A9) 



The reason for the latter simplifications is the following. If 
in = 2 or ui = 6 and d = d* = 4 + m/2 (upper critical di- 
mension), then Bessel functions K v with v = ±| are encoun- 
tered in the integral (Al), which are simple exponentials. c3 
This entails that the required sin gle int e gration s can be done 
analytically to obtain the results ( i49a| )-( 4.50c| ) for the O (e 2 ) 
coefficients. 

For the remaining values of m, i.e., for m = 1,3, 4, 5, the 
required integrals did not simplify to a degree that we were 
able to compute them analytically. However, proceeding as 
explained in Appendix [|, they can be computed numerically. 
In the special cases m = 2 and m = 6, the results of our 
numerical integrations are in complete conformity with the 
analytical ones. 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF THE 
SCALING FUNCTIONS AND E(v) 



Note that for m = 2 the asymptotic form (Bl) is consistent 
with the simple exponential form ( ^\6| ) since $00(2, 5) = 0. 



However, for other values of m, the coefficient (B4) does not 
vanish. For example, $00(6,7) = 1/(8 7r 3 ), in conformity 



with expression (A8) for the scaling function $(u; 6, 7). 



In order to obtain precise results for the integrals 
Jo^im, d*), Ji : z(m,d*), and I\(m,d*), on which the co- 
efficients b<f,(m), b a {rn), and b p (m) depend, we proceeded 
as follows. We split the required integrals as J °° . . . dv — 
Jq° . . . dv + . . . dv, choosing v ~ 9.3. In the integrals 

. . . dv, we replaced the integrands by their asymptotic 
forms obtained upon substitution of $ and/or 5 by their large- 
v approximations $( as ) and s( as ) given in (Bl) and (B2), 
respectively, and then computed these integrals analytically. 
The integrals f^ . . . dv were computed numerically, using 
Mathematicao We checked that reasonable changes of vq 
have negligible effects on the results. The procedure yields 
very accurate numerical values of the requested integrals. The 
reader may convince himself of the precision by comparing 
the so-determined numerical values of the integrals for m = 2 
and m = 6 with the analytically known exact values. 



According to (A4), the scaling function $(u;m, d*) is 
a difference of a hypergeometric function iF 2 and a prod- 
uct of a Bessel function I m /4 times a power. If one asks 
MathematicacJ to numerically evaluate expression (A4) for 
$(u) without taking any precautionary measures, the result 
becomes inaccurate whenever v becomes sufficiently large. 
We found that such a direct, naive numerical evaluation fails 
for values of v exceeding vq ~ 9.5. This is because both 
functions of this difference increase exponentially as v — > 00. 

To cope with this problem, we determined the asymptotic 
behavior of the scaling functions $(u; m, d*) and S(v; m, d*) 
for v — > 00. From the integral representations (3.4) and ( 3.20 ) 
of these functions one easily derives the limiting forms 



<fr(w;m, d) w $( as )(u;m, d) = v 



$oo(m,d) (Bl) 



and 



E{v;m,d) w ■'"»(«; m, d) = v' 2 ^ ^ (m ' % , Hi!) 



I(l-e) 



with 



9 II «i 



e^ii^ll) 
1\\ +ll 



2 2 ( d -™)- 5 7r i ^r(d-2-^) 



(B3) 



r(2=ii3a) 

At the upper critical dimension, the latter coefficient becomes 



$oo(m, d*) 



(B4) 
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Both Sak and GrestH as well as Mergulhao and CarneiroLJ ex- 
plicitly employed the simple form (46) of the sealing function 
in their computations. However, the latter authorsEj did not take 



roil 



advantage of working with the simple function (A8) in the case 
m = 6, performing complicated commutations in the momentum 
representation instead. Sak and GresuJ, on the other hand, did not 
present any details of their calculation for the case m = 6. 
' Mathematica, version 3.0, a product of Wolfram Research. 
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